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The renowned general epidemic process describes the stochastic evolution of a population of indi- 
viduals which are either susceptible, infected or dead. A second order phase transition belonging to 
the universality class of dynamic isotropic percolation lies between endemic or pandemic behavior of 
the process. We generalize the general epidemic process by introducing a fourth kind of individuals, 
viz. individuals which are weakened by the process but not yet infected. This sensibilization gives 
rise to a mechanism that introduces a global instability in the spreading of the process and therefore 
opens the possibility of a discontinuous transition in addition to the usual continuous percolation 
transition. The tricritical point separating the lines of first and second order transitions constitutes 
a new universality class, namely the universality class of tricritical dynamic isotropic percolation. 
Using renormalized field theory we work out a detailed scaling description of this universality class. 
We calculate the scaling exponents in an e-expansion below the upper critical dimension dc = 5 for 
various observables describing tricritical percolation clusters and their spreading properties. In a 
remarkable contrast to the usual percolation transition, the exponents /3 and governing the two 
order parameters, viz. the mean density and the percolation probability, turn out to be different at 
the tricritical point. In addition to the scaling exponents we calculate for all our static and dynamic 
observables logarithmic corrections to the mean-field scaling behavior at dc — 5. 

PACS numbers: 64.60.Ak, 05.40.-a, 64.60.Ht, 64.60.Kw 



I. INTRODUCTION 

The formation and the properties of random structures 
have been a major theme in statistical physics for many 
years. In case the formation of such structures obeys lo- 
cal rules, these processes can often be expressed in the 
language of population growth. It is well known that two 
particular growth processes lead to random structures 
with the properties of percolation clusters. The so called 
simple epidemicprocess (SEP) leads to directed percola- 
tion (DP) dllllli. The SEP is also known as epidemic 
with recovery, Gribov process MM, as the stochastic ver- 
sion of Schlogls first reaction [3, la] or in elementary par- 
ticle physics as Reggeon field theory [gl llol lTli . The so 
called general epidemic process (GEP) |l2l Il3l IT^. also 
known as epidemic with removal, generates isotropic per- 
colation (IP) clusters 15, 16, Ji J^S] and models there- 
fore the universality class of dynamic isotropic percola- 
tion (dip). 

Epidemic models like the SEP and the GEP are rele- 
vant for a wide range of systems in physics, chemistry, bi- 
ology, and sociology. Undoubtedly, the potential of such 
simple models has its limitations because they rely on 
strong simplifying assumptions such as the homogene- 
ity of the substrate To^, isotropy of the infections, im- 
mobility of individuals and so on. However, the transi- 
tion between population survival and extinction of these 
processes is a nonequilibrium continuous (second order) 
phase transition phenomenon and hence is character- 
ized by universal scaling laws which are shared by en- 
tire classes of systems. Near these transitions, simplistic 



epidemic models like the SEP and the GEP are of great 
value, because they are powerful workhorses to study the 
mutual properties of their entire universality class which 
also should include more realistic models. 

The universal properties of DP and dIP are well 
known today thanks to numerous studies of the SEP 
and GEP, respectively. Relatively little is known, on 
the other hand, whether and under what modifications 
these stochastic growth processes allow for first order 
phase transitions between their endemic and pandemic 
states, and by the same token, for tricritical behavior at 
the phase-space boundary between first and second order 
transition. In the context of DP, these questions where 
addressed to some extend by Ohtsuki and Keyes [20| . 
In this paper we will study this intriguing topic in the 
context of dIP by generalizing the famous GEP. 

The standard GEP, assumed for simplicity to take 
place on a lattice can be described with help of the 
reaction scheme 

S'(n)+X(m) X{n)+X{tn), (1.1a) 

X{n) E{n), (1.1b) 

with reaction rates k and A. S, X, and E respectively de- 
note susceptible, ill, and dead (or immune) individuals on 
nearest neighbor sites n and m. A susceptible individual 
may be infected by an ill neighbor with probability k [re- 
action (|l.la|) ]. By this mechanism the disease (henceforth 
also called the agent) spreads diffusively. Ill individuals 
die with a probability A [reaction l)l.lb|l ]. There is no 
healing of infected individuals and no spontaneous infec- 
tion. In a finite system the manifold of states without 
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any infected individual is inevitably absorbing. Whether 
a single infected site leads to an everlasting epidemic in 
an infinite system depends on the ratio k/A. With A 
fixed there is a certain value k = Kc that for all k > Kc 
an eternal epidemic (a pandemic) occurs. The probabil- 
ity P{k) for the occurrence of a pandemic as a function 
of K goes to zero continuously at the critical point Kc- 
The behavior of the model near this critical point is in 
the universality class of dIP. 

Of course, one can conceive many alterations of the 
GEP. Some modifications will lead to models which still 
belong to the dIP universality class whereas other alter- 
ations will produce models belonging to other universal- 
ity classes. As an example for the latter, one might think 
of mobile susceptible or immune individuals which, say, 
diffuse through space. The resulting models are reaction- 
diffusion type models which have nothing to do with 
dip 01 . We are not interested in this kind of alter- 
ations. Rather, we are interested in modifications of the 
GEP that preserve the dIP universality class and which, 
for certain parameter values, allow for tricriticality and 
first order percolation transitions. To be more specific, 
we are interested in the most relevant mechanism lead- 
ing to tricritical and first order dIP. In spirit our model 
is closely related to the canonical model for tricriticality 
in eqilibrium systems, viz. the (j)'^ model with an addi- 
tional (j)^ term where the free energy density is of the 
form / = Tcf)'^ + gicf)"^ + g^cj)^ If 9i is positive then 

higher order terms including the 4>^ term are irrelevant 
and one has in mean field theory a second order tran- 
sition when T passes through zero. Otherwise, however, 
one has an instability and higher order terms are required 
for stabilization of which the 0^ term is the most relevant 
one. For 94 < one has a first order transition at a crit- 
ical value of T that depends on (^4 and f^nd the point 

= .94 — constitutes a tricritical point. Our model to 
be defined in the next paragraph is so that it introduces 
a similar instability in the GEP and via this instability it 
allows for tricriticality and first order percolation. In its 
field theoretic description our model features, compared 
to the standard GEP, an additional higher order term 
which becomes the most relevant stabilizing term when 
the usual GEP coupling vanishes. 

Our modification of the GEP can be defined in sim- 
ple terms. The basic idea is to enrich the reaction 
scheme by introducing weak individuals Y . Instead 
of being infected right away by an ill neighbor, any sus- 
ceptible individual may be weakened with a reaction rate 
/i by such an encounter. When the disease passes by 
anew, a weakened individual is more sensitive and gets 
sick with a rate v > k. In the following we refer to this 
model as the generalized GEP (GGEP). In addition to 
the reactions the GGEP is described by the reac- 

tions 

S{n)+X{m) ^ y(n) + X(m), (1.2a) 
y(n)-fX(m) X{n) + X{ra). (1.2b) 

As we go along, we will show that the occurrence of the 



weak individuals gives rise to an instability that can lead 
to a discontinuous transition and compact (Eden |2^i2^]) 
growth of the epidemic if v is greater than a critical value 
i'c(K,/i). In the enlarged three-dimensional phase space 
spanned by k, /i, and u with fixed A, there exists a critical 
surface associated with the usual continuous percolation 
transition and a surface of first order transitions char- 
acterized by a finite jump in the probability P(k, /x, v) 
for the occurrence of a pandemic. These two surfaces of 
phase transitions meet at a line of tricritical points. 

The focus of this paper lies on the universal proper- 
ties of the GGEP near this tricritical line. By using the 
methods of renormalized field theory we work out a scal- 
ing description of the new universality class of tricritical 
dynamic isotropic percolation (TdIP) to which the tri- 
critical GGEP belongs. We study a multitude of static 
and dynamic observables that play important roles in 
percolation theory. In particular we calculate the crit- 
ical exponents describing the scaling behavior of these 
observables below 5 dimensions as well as logarithmic 
corrections to the mean-field scaling behavior in 5 di- 
mensions. 

The outline of our paper is as follows. In Sec. ITU we 
consider the GGEP in a mean-field theory. As the main 
result of Sec. ^ the mean-field analysis will reveal the 
structure of the phase diagram. With the aim of study- 
ing the effects of fluctuations, we condense the principles 
defining TdIP into a field theoretic model in Sec. IIIII 
In Sec. IIVI we work out the scaling properties of static 
aspects of TdIP. Section treats the dynamic scaling 
properties. Goncluding remarks are provided in Sec. IVII 
There is one Appendix in which we sketch the calcula- 
tion of a parameter integral that is helpful in computing 
Feynman diagrams. 

II. MEAN-FIELD THEORY 

A mean-field description of the GGEP can be formu- 
lated by treating the reaction equations (|l.ll) and (|1.2|l as 
deterministic equations without fiuctuations. This deter- 
ministic approximation leads to the system of differential 
equations 

nn(n) 

S{n,t) = -{k + ii)S{-a.,t) X{m,t), (2.3a) 

m 

nn{n) 

Y{n,t)=(iiS{n,t)-vY{n,t)) ^X(m,i), (2.3b) 

m 

Tin(n) 

i:(n,i) = (^K5(n,t)-f z^r(n,t)) ^ X{m,t) 

m 

-AA(n,t), (2.3c) 
E{n,t) = \X{n,t) (2.3d) 

governing the dynamics of the different kinds of individ- 
uals. Here, X)^*""' denotes summation over the nearest 
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neighbors of n. At each lattice site there is the additional 
constraint 



S + X + Y + E = 1 



(2.4) 



Thus, S, X, Y, and E can be interpreted as the prob- 
ability of finding the corresponding state at a site n. 
Note that the processes can only proceed at places where 
the probability to find ill individuals in the neighbor- 
hood is not zero. We use the canonical initial condition 
that all sites of the initial state are susceptible except 
for the site at the origin which is ill, Xo(n) — 5nO and 
Yo{n)=Eo{n) = 0. 

Equations (|2.3a|l and l|2.3bl) are are readily integrated. 
We obtain 

/ nn(n) , 

S{n, t) = 5o(n) exp - p ^ E{m, t) ) (2.5) 



and 

/x5(n, t) + {p- y)Y{n, t) = [pSo{n) + {p - y)Yo{n)) 

nn(n) 



exp I —u E{m.^t) 



(2.6) 

where we have defined p = k + p. The time scale A has 



been set to unity for simplicity. Equation l)2.3d|l together 
with the constraint H2.4|l leads finally to the mean-field 
equation of motion of the GGEP, 



Ein,t) 



1 - E{n,t) 
P 



K 



p-v 



exp 



p-v 



■ exp 



nn(n) 

V ^ E{m,t) 



nn(n) 

p J2 E{,n,t))\So{n). (2.7) 



In the asymptotic regime \n\,t — > oo one can neglect 
time and space dependence in (|2.7|l and use the approxi- 
mation X)^''"'' E!{m., t) w zE{n, t), z being the coordina- 
tion number of the lattice. Hence the asymptotic values 
of E are the solutions of the equation of state 



E = 



1 



e 



-zpE\ _ 



= --f{E). 



p — V p — V 

(2.8) 

By setting p = corresponding to = 0, one obtains 
the equation of state of the usual GEP 



E = 1 



(2.9) 



with K, — — 1/z determining the second order phase 
transition corresponding to ordinary isotropic percola- 
tion. For K < Kc, Eq. (|2.9|l has only the solution E = 
which means that the disease does not percolate, i.e.. 



is endemic. In the other case, k > Kc, a stable solu- 
tion E > arises, signalling the percolative pandemic 
character of the disease. These types of solutions exist 
also in the full equation of state H2.8|l . Using k < p, 
one demonstrates easily that f{E) increases monotoni- 
cally from zero to one in the interval < £' < oo. Thus, 
i? = is always a solution of Eq. H2.8|) . It follows from the 
equation of motion (|2.7|) that only solutions of Eq. H2.8|) 
with f'{E) < 1 are stable. Because /'(O) = zk, a stable 
percolating solution E > exists always for k > Kc- The 
existence of more than one nontrivial solution requires 
necessarily that f"{E) = at least for one value E > {). 
From Eq. H2.8|l we derive the inequality 

f"{E)/z^ < {^pv - k{p + v)^ exp{-zEmax{p, v)) . 

(2.10) 

With f"{E) < Ofori;> 1 and f"{0)/z^ ^ pv-K{p + v) 
we therefore find 

^ (2.11) 



1 1 
- > - 

K p 



as the necessary and sufficient conditions for the existence 
of a second locally stable solution E > Q besides E = Q 
with a discontinuous transition between the two. The line 
of tricritical points where the first and the second order 
transitions meet in the K-p-v phase space is determined 
by 1/k = 1/ p + 1/v = z. 

In the following we focus on the phenomena arising 
near the tricritical line. In this region of the phase space, 
E(n,t) is, except for a microscopic region around the ori- 
gin n = of the position space, small and slowly varying. 
Hence, we may approximate Eq. I|2.7|l by the determin- 
istic reaction-diffusion equation 



£;(x,t) = DV^Eix.t) - A 



TE{yi, t) 



f 



E{^,tf + ^-E{^,t) 
6 



(2.12) 



where D = 



with a being the lattice constant, t = 



1 — k/kc, f — Kpv{l/ p+ 1/v— 1/k)/k^, and g = {pv{p + 
V + k) — k{p + i^)^)/k|? w pv/n^. As a consequence of 
Eq. Ij2.3d|l we have the constraint E{x,t) > 0. 

Holding the positive coupling g constant, Eq. (|2.12|) 
comprises only two tunable parameters, viz. r and / 
and hence the dimensionality of the phase diagram is 
reduced to two. As follows from the different types of 
solutions of Eq. H2.12|l to be discussed in a moment, this 
two-dimensional phase diagram features a line of second 
order transitions (A- line), a line of first order transitions 
and a tricritical point determined by r = / = separat- 
ing the two lines of transitions. In addition, there are 2 
spinodal lines. The entire phase diagram is depicted in 
Fig. [2 

Equation (|2.12|l has besides the trivial solution E = 
0, which is stable for r > 0, a nontrivial locally stable 
stationary homogeneous solution E = A with 



A= (^9p/4-6gT~3f/2 



(2.13) 
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FIG. 1: The mean-field phase diagram. 



which is physical only if A > 0. For / > one has 
therefore a continuous transition from E — {t > 0) to 
E ^ A !=i 2\t\ /f {t < 0). However, for / < there 
is a certain value t — between the spinodals r = 
3/^/8.9 and r = for the discontinuous transition from 
E = to E = A. The values of may be determined 
by studying a travelling wave solution of the equation 
of motion (|2.12l) that describes the infection front of a 
big expanding sherical cluster of the epidemic. Such a 
solution is given by 



E{x, t) = a(i - ta.Bh{b{x - vt))^ /S 
X{x,t) = E{x,t) 



vbA 



2[cosh(fo(a; - vt))]' 



(2.14a) 
(2.14b) 



with 



"exact" approach which, as a first step, consists of re- 
formulating the microscopic master-equation for the re- 
actions (11.111 and p.2|l as a bosonic field theory on the 
lattice. The next and pivotal step in this approach is 
to take the continuum limit. Albeit the "exact" ap- 
proach with a naive continuum limit leads to a consis- 
tent dynamic functional for the GGEP (after deleting 
several irrelevant terms) this approach must be cautioned 
against. Strictly speaking, one has to use Wilson's statis- 
tical continuum limit |2q in the renormalization theory 
of critical phenomena. This procedure consists of suc- 
cessive coarse graining of the mesoscopic slow variables 
(order parameters and conserved quantities as functions 
of microscopic degrees of freedom), the elimination of 
fluctuating residual microscopic degrees of freedom, and 
a rescaling of space and time. In general, microscopic 
variables do not qualify as order parameters. Alarm- 
ing examples are the pair contact processes (PGP and 
PGPD), where a naive continuum limit of the micro- 
scopic master-equations leads to untenable critical mod- 
els. Therefore, we devise our field theoretic model repre- 
senting the TdIP universality class using a purely meso- 
scopic stochastic formulation based on the correct order 
parameters identified through physical insight in the na- 
ture of the critical phenomenon. Hence, our dynamic 
response functional stays in full analogy to the Landau- 
Ginzburg- Wilson functional and provides a reliable start- 
ing point of the field theoretic method. 



A. Langevin equation 



A5 



A. 



I3XD 



(x/9/V4-6ffT + //2 



(2.15) 

The condition X(x,t) > requires v > 0. The first order 
transition from E = to E = A a,t is therefore defined 
by the phase equilibrium condition v = leading to 

rd = ^. (2.16) 
3.9 

E jumps at T — from zero to the value 2\f\ /g. 



III. FIELD THEORETIC MODEL 

In this section we will derive a dynamic response func- 
tional J5, 26, 27] for the GGEP based on very general 
arguments alluding to the universal properties of TdlP. 
First we will distill the basic principles of percolation 
processes allowing for tricritical behavior. Next, we cast 
these principles into the form of a Langevin equation. 
Then we refine the Langevin equation into a minimal 
field theoretic model. 

As an alternative avenue to a field theoretic model for 
the GGEP one migtht be tempted to use the so-called 



The essence of isotropic percolation processes can be 
summarized by four statements describing the universal 
features of the evolution of such processes on a homoge- 
neous substrate. Denoting the density of the agents (the 
infected individuals) by n(r, t) and the density of the de- 
bris (the immune or dead individuals) which is propor- 
tional to the density of the weakened substrate by m(r, t), 
these four statements read: 

(i) There is a manifold of absorbing states with n = Q 
and corresponding distributions of m depending on 
the history of n. These states are equivalent to the 
extinction of the epidemic. 

(ii) The substrate becomes activated (infected) de- 
pending on the density of the agents and the den- 
sity of the debris. This mechanism introduces 
memory into the process. The debris ultimately 
stops the disease locally. However, it is possi- 
ble that the activation is strengthened by the de- 
bris through some mechanism (sensibilization of the 
substrate). 

(iii) The process (the disease) spreads out diffusively. 
The agents (the activated substrate) become deac- 
tivated (converted into debris) after a short time. 
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(iv) There are no other slow variables. Microscopic de- 
grees of freedom can be summarized into a local 
noise or Langevin force C(r, t) respecting the first 
statement (i.e., the noise cannot generate agents). 

The general form of a Langevin equation resembling 
these statements is given by 

X-^n ^V^n + R{n,m)n + ( , (3.1a) 



m(r,t) = A / n{r,t')dt' , 



(3.1b) 



where A is a kinetic coefficient and the Gaussian noise 
correlation reads 



C(r,t)C(r',i') = X-^Q{n, m) n{r, t) 5{y - r') 5{t - t') 

- \-^an{r, t) V2,5(r - r') 6{t - t') 

-I- Q'{n, m)n(r. t)n(r, t') S(r - r') -\ . 

(3.2) 

The first row in Eq. H3.2|l represents time-local reaction 
noise. The second row describes noise originating from 
diffusion and the last row shows an example of possi- 
ble time-non-local noise (quenched noise) that may be 
acquired through random disorder or through the elimi- 
nation of microscopic slow variables, e.g., fluctuations of 
the debris. The structure of the 3 terms is so that they 
respect the absorbing state condition. Of course many 
further contributions to Eq. (|3.2|) are conceivable (hence 
the ellipsis in the third row) including non-Markovian 
and also non-Gaussian noise. We will see below, that all 
these terms turn out being irrelevant, and that only the 
simplest form of reaction noise contributes the minimal 
field theoretic model. 

The dependence of the rate R{n, m) on the density 
of the debris m(r, t) describes memory of the process 
mentioned above. We are interested primarily in the be- 
havior of the process close to the tricritical point, where 
n and m are small allowing for polynomial expansions 
i?(n, m) = —T—an — fm—gm^/2 + . . ., Q(n, m) = 7+- • • , 
and Q' {n, m) = 7' -I- • • • . The justification for the trunca- 
tion of the expansions will be given later by IR relevance- 
irrelevance arguments. As long as / > 0, the second 
order term gm? of the rate R is irrelevant near the tran- 
sition point and the process models ordinary dynamic 
isotropic percolation 0| . We permit both signs of 
/ so that our model accounts for sensibilization (weak- 
ened substrate) and allows for compact (Eden) spreading. 
Consequently we need the second order term for stabi- 
lization purposes, i.e., to limit the density n to finite 
values. 



B. Dynamic response functional 

In order to apply the renormalization group (RG) and 
field-theoretic methods H^H^I, it is convenient to use the 
path-integral representation of the underlying stochastic 



process n{r,t) |2^|2^|23. With the imaginary- valued 
response field denoted by n{r,t), the generating func- 
tional of the Green functions (connected response and 
correlation functions) takes the form 

W[H,H] ^\nj V[n,n] exp|-J^[H,n] 

+ J d'^x J dt{Hn + Hn)Y (3.3) 

The generating field H{r, t) corresponds to an addi- 
tive source term for the agent in the equation of mo- 
tion H3.1|) . Therefore, the response function defined by 
a functional derivative with respect to H{r,t) describes 
the influence of a seed of the agent at {r,t). The dy- 
namic functional [n, n] and the functional measure 
r'[n,n] cx Yl^^[dn{r,t)dn{r,t)/2TTi) are defined using 
a prepoint (Ito) discretization with respect to time [2^ . 
The prepoint discretization leads to the causality rule 
9{t < 0) = in response functions. This rule will 
play an important role in our diagrammatic perturba- 
tion calculation because it forbids response propagator 
loops (see below). Note that the path integrals are 
always calculated with the initial and final conditions 
n(r, —00) = n{r, 00) = 0. 

The stochastic process defined by Eqs. (|3.1|) and 13. 2|) 
leads via the expansions of i?, Q, and Q' to the prelimi- 
nary dynamic response functional 

J' ^ J d'^x |a j dt n{\-^dt -W^ + T + an + fm 



9 2\ 1-2 /'Y7~^2 

—m ]n nn — an Vn 



' r r _ 1 2 
\ dt nn 



(3.4) 



As we have remarked above, the term proportional to g 
may be neglected only if the coupling / is positive defi- 
nite. In this case Eq. (|3.4|l reduces to the response func- 
tional of usual dynamic percolation 0, As in all 
models with an absorbing state transition, the functional 
J^' includes a redundant variable which has to be removed 
before any application of relevance- irrelevance arguments 
since it has no definite scaling dimension. This redundant 
variable is connected with the rescaling transformation 



n— n— >6 ^n, a— 



b-' 



9, 7 



(3.5a) 

67, 
(3.5b) 



which leaves invariant. Exploiting this invariance, we 
we may set 7=1 which fixes the redundancy. Of course, 
this is justified only if 7 is a finite positive quantity in 
the region of interest of the phase diagram. 

For the steps to follow, we need to know the naive di- 
mensions of the constituents of J^' after the removal of 
the redundant variable. As usual, we introduce a con- 
venient external length scale so that r ~ fi~^ and 



Xt ^ fj, ^. Exploiting that J'' has to be dimensionless, 
we readily find 

nr^ fi^ , fj.'^-^ , TO - ^'^-^ , (3.6a) 

r~/i2, <?^m'<5-'*', (3.6b) 

a~/-^ 7'-/-''. (3.6c) 

The highest dimension d at which any of the finite and 
positive couplings becomes marginal (acquires a vanish- 
ing naive dimension) corresponds to the upper dimension 
dc of the theory. This dimension separates trivial mean- 
field critical behavior for d > dc from non trivial behavior 
in the regime where d < dc and where the relevant vari- 
able T is small. Thus, if / is finite and positive it follows 
that dc — 6. Then g, a, 7', and a have negative naive 
dimensions and are therefore irrelevant. The correspond- 
ing terms in the response functional H3.4|l vanish at the 
critical fixed point, and the response functional displays 
the asymptotic symmetry f^^/^n{r,t) ^ ~ f^/^Tn{r, ~t) 
. The resulting response fmictional is that of the GEP 

However, if / is zero as it is at the tricritical point, 
we must use g to fix the upper critical dimension which 
leads to = 5. The dimensions of a, 7', and a are nega- 
tive near dc- Thus, diffusional and quenched randomness 
of the noise is irrelevant here as they are for the GEP. 
Dimensional analysis also justifies the truncation of the 
expansions of P, Q, and Q' , as well as the elimination of 
other terms. All higher order terms are irrelevant in the 
IR renormalization group approach because they carry- 
ing a negative naive dimension near d — dc- Collecting, 
we obtain the dynamic response functional 

J = J d'^xX j dth(^X-^dt-V^ + T 

+ fm+^m^^^fi^n (3.7) 

as our minimal field theoretic model for the TdIP uni- 
versality class. We shall see as we move along, indicating 
the consistency of our model, that the proper elimination 
of IR irrelevant terms has led to an UV renormalizable 
theory at and below five dimensions. 

Note that, contrary to dIP, the functional (|3.7|) does 
not have an asymptotic symmetry that relates the two 
fields h and n. Thus, these fields will get independent 
different anomalous dimensions leading to two different 
order parameter exponents. 

Before we go on, we extract the diagrammatic elements 
implicit in J^. These elements will play a central role later 
on when we calculate the Green functions perturbatively. 
As usual, the Green functions are the cumulants of the 
fields n and n which correspond in graphical perturbation 
expansions to the sums of connected diagrams. Their 
actual calculations are performed most economically in 
a time-momentum representation. To this end we will 
use the spatial Fourier transforms of the fields n and n 
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^ ^ 



FIG. 2: The diagrammatic elements implied in the dynamic 
functional J. 

defined via 

n(r,t)= / e^'i '-nq(t) , n{r,t) ^ ( e'^ '^n^it) , (3.8) 

where /q • ■ ■ := (27r)~'' / d'^q ■ ■ ■ . In the time-momentum 
representation we can simply read off the diagrammatic 
elements from the dynamic functional. The harmonic 
part of comprises the Gaussian propagator 

(«q(<)nq' it'))o = (27r)'^G(q, t - t') (5(q + q') , (3.9a) 
G(q, t) = 6{t - t') exp ( - A(r + c^)t^ , (3.9b) 

where (■ • • )o indicates averaging with respect to the har- 
monic part of the dynamic functional (|3.4() . The non- 
harmonic terms give rise to the vertices A, —X^f9{t — t') 
and —X^g9{t — t')9(t — t"). AU four diagrammatic ele- 
ments of J are depicted in Fig. [5] 

C. Quasi-static model 

As mentioned earlier, we are interested in the dynamic 
as well as the static, i.e. t — > 00, properties of tricriti- 
cal percolation. Of course, we could base our entire RG 
analysis on the full dynamic functional as given in 
Eq. (|3.4|l . Then we could extract the static behavior 
from the dynamic behavior in the end by letting t — > 00. 
This would mean, however, that we had to determine 
all the required renormalizations from dynamic Feyn- 
man diagrams composed of the diagrammatic elements 
listed in Fig. |2] Fortunately, there is a much more eco- 
nomic approach possible here which is based on taking 
the so-called quasi-static limit. We will see shortly, that 
the perturbation theory simplifies tremendously in this 
limit. All but one renormalization factors can be calcu- 
lated using this much simpler approach. Only for the 
one remaining renormalization we have to resort to the 
dynamic response functional . Taking the quasi-static 
limit amounts to switching the fundamental field variable 
from the density of the agents to the density of the de- 
bris TO(r) := TO(r, 00) — X J^^dtn{r,t) left behind by 
the epidemic. Then, the static properties of TdIP can be 
studied via the Green functions of the debris density. In 
particular the response functions (J^^ TO(ri)H(0, 0)) and 
their connected counterparts will be important for our 
analysis because they encode the static properties of the 
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-A = G(q) 




-f <- 



FIG. 3: The diagrammatic elements implied in the quasi- 
static Hamiltonian Ti.. 

percolation cluster of the debris emanating from a seed 
localized at the origin at time zero. 

After this prelude we now formally take the quasi- 
static limit of the dynamic functional. The structure of 

is so that we can directly let 



n(r, t) s{r) , m(r) = A / dt n{r, t) — > s(r) . 

J —oo 

(3.10) 

This procedure leads us from to the quasi-static Hamil- 
tonian 



fs+^-s^-\i)s. (3.11) 



It is easy to see that Ti generates each Feynman dia- 
gram that contributes to (flj m(rj) n(f ^ , 0)). Stand- 
ing alone, however, this Hamiltonian is not sufficient to 
describe the static properties of tricritical isotropic per- 
colation (TIP). As a remainder of its dynamical origin Ti, 
must be supplemented with the causality rule that for- 
bids closed propagator loops. 

The propagator of the quasi-static theory follows from 
Eq. H3.11|l [or likewise from Eq. (|3.9|) ] as 



SqSqOo = (2^)''G(q) <5(q + q') , G(q) = 



1 



r + q" 

(3.12) 

As far as vertices are concerned, TL implies the three- 
leg vertices 1 and — / and the four- leg- vertex —g. 
The quasi-static propagator and vertices are shown in 
Fig. 13 Knowing all the diagrammatic elements, one 
can straightforwardly check in explicit graphical pertur- 
bation expansions that the quasi-static Green functions 
calculated with Ti are equal diagram for diagram to the 
zero- frequency limit of the dynamic Green functions (the 
Green functions of the time-integrals of n(r, t)) calculated 
with J. 

Before embarking on our RG analysis, we finally men- 
tion the naive dimensions of the quasi-static fields. These 
are given by 



d-A 



(3.13) 



IV. STATIC SCALING PROPERTIES 

Now we will study the static properties of TdIP, that 
is the properties of TIP. The dynamic properties of TdIP 
will be addressed later on. 



A. Diagrammatics 

In principle, we could extract the properties of TIP 
directly from the correlation functions 



(s(ri) • • • s(rAr) s{vn+i) ■ ■ ■ s{vj^,j^)) 



I?[s,3]s(ri)...s(r^^^)e-«P-1 



(4.1) 



However, since our model is translationally invariant, it is 
much more convenient to use the vertex functions ^ ^ 
instead which are related to the connected counterparts 



Gat.^IW) =(.<ri)...s(rjv) 



(conn) 



(4.2) 



of the correlation functions via Legendre transformation 
of their generating functional js^l . Graphically, the ver- 
tex functions consist of amputated one-line irreducible 
diagrams. As usual in determining the renormalizations, 
we can restrict ourselves to the superficially divergent 
vertex functions, i.e., those vertex functions that have a 
non-negative //-dimension at the upper critical dimension 
dc- A simple dimensional analysis shows that only the 
vertex functions that correspond to the different terms 
of the Hamiltonian (|3.11l) . viz. Fu, Fi ^, r2,i, and Fi,3, 
are primitively divergent. Thus, the theory is renormal- 
izable by additive and multiplicative renormalization of 
the fields and the parameters of the theory. 

Throughout, we will use dimensional rcgularization 
to calculate the Feynman diagrams constituting the re- 
quired vertex functions, i.e., we will compute the dia- 
grams in dimensions where they are finite for large mo- 
menta and then continue the dimension analytically to- 
wards dc- This procedure converts the logarithmic large- 
momentum singularities of a cutoff rcgularization into 
poles in the deviation e = dc — d from dc- However, poly- 
nomial large-momentum singularities that require addi- 
tive renormalizations are unaccounted for in dimensional 
rcgularization. In order to discuss such additive renor- 
malizations, we will occasionally use a cutoff rcgulariza- 
tion with a large momentum cutoff A. 

At this place it is worth stressing that real critical 
systems not involve any UV-divergencies because all in- 
verse wavelengths of fluctuations have a physical cutoff 
A. On the other hand, critical systems suffer from IR- 
divergencies. However, if and only if we use a reliable 
field theory, we can formally transform IR-divergencies 
into UV-divergencies by a simple rescaling. This way we 
can learn about IR-scaling properties of the critical sys- 
tem indirectly via the UV- renormalizations. A correct 
and reliable statistical field theory constitutes what Wil- 
son (25 calls a logarithmic theory free of length scales. 
Only within such a theory it makes sense to apply the 
techniques of renormalized field theory to critical sys- 
tems. 





(a) 



(b) 



FIG. 4: One-loop diagrams contributing to the renormaliza- 
tion of Fi,! and Fi,2 if a momentum-cutoff regularization is 
used. In dimensional regularization, these diagrams are finite. 



1. Divergent One-Loop Diagrams 

Two divergent one-loop diagrams can be assembled 
from the diagrammatic elements listed in Fig. O These 
diagrams, which contribute to the vertex functions Fi^i 
and Fi^2, are shown in Fig.^ It is easy to see that they 
are linearly divergent for d = dc = 5 (i.e., the /i dimen- 
sion is 1). No logarithmic divergencies arise at one- loop 
order. Thus, the theory can be renormalizcd to one-loop 
order by additive renormalizations 

f^f = ^i'{v-bAfi'''g), (4.3a) 

r ^ f = r - aAfi^'^f = (t - aAv) + abk^^x^'^^g , 

(4.3b) 

where the open circles indicate unrenormalized quantities 
and where a and b are positive constants. 

Next we will switch to dimensional regularization for 
convenience. In this method all polynomial divergencies 
arising from the large cutoff are formally set to zero, 
and hence the two diagrams in Fig. 0] become finite at 
dc- Thus, the additive renormalizations Eq. (|4.3I) be- 
come formally superfluous in dimensional regularization. 
However, it has to be emphasized that dimensional regu- 
larization is only a formal trick. Physically the additive 
renormalizations are always present and we need the in- 
teraction term proportional to the coupling constant f t o 
renormalize the theory contrary to the claim of Ref. |2Clj| . 
Using dimensional regularization we have to keep in mind 
that these additive renormalizations do exist and that the 
critical "temperature" is shifted by a term linear in the 
crossover variable u: fc = aKv — abK'^^.~'^'^g-\-0{g^). This 
fc is formally set to zero in dimensional regularization. 



2. Two-Loop calculation 

Since there are no e poles at one-loop order, we have to 
proceed to higher orders in perturbation theory to find 
non-trivial critical exponents. We will see that e poles do 
occur at two-loop order and that the two-loop diagrams 
will lead us to anomalous contributions to the critical 
exponents of order e. 

We start with the self-energy. The two-loop diagrams 
contributing to the renormalization of Fi.i are listed in 
Fig. 13 In the following we will use a compact notation 
for the diagrams. For example, jSJi) refers to diagram 




<2> 

(b) 




(c) 



FIG. 5: Two-loop diagrams contributing to the renormaliza- 
tion of Fi.i. 



(a) of Fig. [S] For the momentum integrals occurring in 
the diagrams we will use the abbreviations 



kin 



qi:q2 



(qf + 1) (q| + l) ((qi+q2)2 + l)' 



(4.4) 

This has the benefit that we can write the divergent parts 
of the diagrams in a compact form. For the self-energy 
we have 



0i) 



fr-^(/ii2r-^W), (4.5a) 



(4.5b) 



The Ikim can be calculated very efficient with help of 
the parameter integral 

I{a,b;c) 



(q2 + a) (q2 + b) ((qi + qa)^ + c) 



2 • 



(4.6) 

Using dimensional regularization, we find the e expansion 
result 

l{a,b;c)^-^{a + b-c) + 0{e') (4.7) 

for the parameter integral. This calculation is sketched 
in the Appendix. In Eq. (|4.7I) we used the shorthand 
G, = F(l + e/2)/(47r)''/2 for convenience. By taking 
derivatives with respect to the parameters a, b and c, we 
get the singular parts 



22 



'112 



-2/i 



13 



3e 



(4.8) 



of the required original integrals. Collecting we find 
Ti.i - (r + q2) - dSh) - ®) - (H:) 



1 + ^,-) 2 
lOe 



(4.9) 



for the singular part of the inverse response function. 
Here, we introduced = for notational conve- 

nience. 

Now we turn to the vertex function Fi^2. Its two-loop 
contributions are shown in Fig. |S1 We obtain 

&) - ®) = e) = Pf) = fgr-'hia , (4.10a) 
int) = 2 • lli) = int) = (Eh) = 2/5r~^/i22 . (4.10b) 




(a) 



(b) 




(c) 



(d) (e) (f) 






(g) 



(h) 



FIG. 6: Two-loop diagrams contributing to the renormaliza- 
tion of ri,2. 





(a) 



FIG. 7: Two-loop diagrams contributing to the renormaliza- 
tion of r2,i. 



Summing up the individual terms we get 



(4.11) 



For the vertex function r2,i we have to calculate the 
diagrams in Fig. With 



we obtain 



13 



-2(7T-^/i22, (4.12) 



i 2,1 = I i T 



(4.13) 



It remains to consider Fi.s. The diagrams in Fig. |H1 
lead to 



2 • (Et) = 2 • ®) = (p) = 3.gV-^/ii3 , 
2 • (|Ht) = B) = 2 • (Elf ) = 6g^T-'Ii22 ■ 



From these expressions we get 

ri,3 



, l8A,g _ 

1 T 



9- 



(4.14a) 



(4.15) 



For completeness, we conclude our quasi-static pertur- 
bation theory by briefly returning to a cutoff regulariza- 
tion. In cutoff regularization there are two additional 
singular two-loop diagrams, see Fig. ^ These diagrams 





(d) 



(e) 



(f) 



FIG. 8: Two-loop diagrams contributing to the renormaliza- 
tion of Fi 3. 





(a) 



(b) 



FIG. 9: Two-loop diagrams that contain insertions of the 
one- loop diagram shown in Fig.jljb). In dimensional regular- 
ization, these two-loop diagrams are finite. 



have divergent insertions of the singular one-loop dia- 
gram (0Jd). Hence, the diagrams of Fig. O are finite in 
dimensional regularization. However, in the more phys- 
ical cutoff regularization they diverge linearly with the 
cutoff. These divergencies are ultimately cancelled by 
the additive rcnormalizations, Eq. 1)4. 3|) . 



B. Renormalization 

Next we absorb the e-poles, into a reparametrization 
of the fields and the parameters of the theory. For the 
quasi-static fields we employ the renormalizations 



(4.16) 



For the parameters of the quasi-static Hamiltonian (|3.11|) 
we use the scheme 



T^T = Z^\ZrT + Yv^), 



where 



Zx = (ZZ) 



1/2 



(4.17a) 
(4.17b) 
(4.17c) 

(4.18) 



It follows from Eq. (|4.1()|) that the vertex functions are 
renormalized by 



N,N 



N,N ~ ^ ^ N,N 



(4.19) 
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Using Eq. (|4.19|l together with our two-loop results H4.9|l . 
gUl), imnj), and we find 

Z=l + 0(u^), ^=1 + + 0(u2), 

he 5e 

(4.20a) 

Z„ = 1 + + 0(^.2) , = 1 + 12!^ + 0(^2) ^ 

(4.20b) 

= 1 + - + 0(ii2) , y = — + 0(u) , (4.20c) 
e Is 



for the renormalization factors. Equation H4.20a|) implies 
Zx = \-^^+0(u^). (4.21) 

C. Renormalization group equation 

In order to explore the scaling properties of tricriti- 
cal percolation we now set up a renormalization group 
equation (RGE). This can be done in a routine fashion 
by exploiting that the bare (unrenormalized) quantities 
must not depend on the arbitrary mesoscopic length scale 
\x~^ introduced in the course of the renormalization. In 
particular the bare Green functions must be independent 
of /i, i.e., 



N,N 



= 0, 



(4.22) 



where denotes /^-derivatives at fixed bare parame- 
ters. Switching from bare to renormalized quantities, the 
identity 14.22|l translates into the RGE 



(4.23) 



Here, I?^ stands for the RG differential operator 

T^ti = + i'^'^r + v'^Kyr)dr + VK^dy + /3„(9„ (4.24) 

that features the Gell-Mann-Low functions 

Pu = fJ.df,\gU^ (^-2e + ^7 + ^7 - w , (4.25a) 

VKy = n df,\gV = ( -e + 7 + ^7 - Ivjv, (4.25b) 



TKr + v'^Kyr = fJ' dfiloT = ( - 7r ) T - Jtvv'^ 



and the Wilson functions 

j = fid^\^\nZ, 7 /i(9^|glnZ, 
7i = M S^loln^i , i^u,v,T. 



(4.25c) 



(4.26a) 
(4.26b) 



From Eq. (|4.25a|l we know that the functions /?„ and 
VKy begin with the zero-loop terms —2eu and —ev re- 
spectively. The higher order terms are determined by 



the Wilson functions. The particular form of these 
functions can be straightforwardly extracted by using 
7... = /i9^|glnZ... = /39„lnZ.... In minimal renor- 
malization the Z-factors have a pure Laurent expan- 
sion with respect to e, i.e., they are of the form Z... — 

1 + Z.^.'Ve + ^-^'Ve' + • • • • Thus, recursively in the loop 
expansion, the Wilson functions also have a pure Laurent 
expansion. Moreover, because the Wilson functions must 
be finite for e — > 0, all the e-poles in this expansion have 
to cancel (this provides a valuable check for higher order 
calculations). Hence, we can obtain the Wilson functions 
readily from the formula 7... = —2uduZ'^^\ The same 
argumentation also leads to ^rv — —2Y^^'> — 2uduY'-^\ 
where Y'^^^ is the coefficient of the first order term in the 
Laurent expansion Y = y /e + --- of the y-factor. Us- 
ing this prescription, we derive from our renormalization 
factors (|4.20|) that the Wilson functions are given by 



7 = — + 0(w^), 

7„ = -36u + 0{u^) , 
7^ = -2u + 0{u^), 

From these results we get 



^^^^ + 0{u^), (4.27a) 
5 

7„ = -20u + O(u2), (4.27b) 

^-3 + 0{u) . (4.27c) 



(3u ^ {-2e + 45u + 0{u^))u, 
Hv = -£H — + Oil/) , 



11m „, 2\ 

+ 0(U ) , K-ar 



3 + 0(u) 

for the Gell-Mann-Low functions (|4.25|) . 



(4.28a) 
(4.28b) 

(4.28c) 



D. Scaling properties 

1. General scaling form 

Next we solve the RGE (|4.23f) by using the method 
of characteristics. The strategy behind this method is 
to introduce a single fiow parameter i that allows to re- 
express the partial differential equation (|4.23|) as an or- 
dinary differential equation in terms of t. This equation 
then describes how the Green functions behave under a 
rescaling 



/i — > /i(^) — 



(4.29) 



of the inverse length scale ji. The characteristic for the 
dimensionless coupling constant u is given by 



' M 



/3„(u(€)), u{l)=u. 



(4.30) 



With the exception of the characteristic for r, the re- 
maining characteristics are all of the same structure, viz. 



dl 



(4.31) 



11 



Here, Q is a placeholder for X, X, and v, respectively, q 
is a placeholder for respectively 7, 7, and k^. The initial 
conditions pertaining to Eq. ^A.[M\ are X{1) — X{1) = 1 
and v{\) = V. The long-length scale behavior of TIP 
corresponds to the limit — > 0. In this limit the RG 
flows to a fixed point determined by the stable value u* 
of u satisfying /3„(u*) = 0. We find that this value is 
given by -u, = 2e/45 + 0{e^). 

For the remainder of Sec. II V 151 we exclusively consider 
dimensions less than 5. We will turn to the case d = 5 in 
Sec. lIVEl In the vicinity of the fixed point the solution 
of the RGB with the characteristics is fairly straightfor- 
ward. We are confronted, however, with the slight com- 
plication that r itself is not a scaling variable as can be 
seen from Eq. (|4.25c|) . In order to diagonalize the flow 
equations in near we switch from t to 



■ V . 



2 ^ 



(4.32) 



where Kt* = k,t{U),), Kv* = and so on. It can 

easily be checked that a is governed by the flow equation 



dln^_ 



(4.33) 



i.e, that a is a true scaling variable. In the dimensions of 
interest here the solutions of the characteristics for the 
scaling variables are of power law form and we obtain 



xG^f^{{v},at\vl''\u^,^i£), (4.34) 



where we have omitted nonuniversal amplitude factors. 
The various exponents appearing in Eq. (|4.34|) are given 
by 



■n 


= liu*) 


00 

= 225^^ 




(4.35a) 


V 


7("*) 


28 

= --^'- 


fO(e2), 


(4.35b) 






22 

)= 225" 


+ 0(e'), 


(4.35c) 


K2 


— 




-0{s'). 


(4.35d) 



Supplementing the solution (|4.34l) with a dimensional 
analysis to account for naive dimensions we obtain the 
scaling form 

G^^^({r},r,.,.,M)=A-M(''-''™ 

X Fj,j^{{efir}, fi-'v/e'^/n , (4.36) 

where 



and where the ^ are, up to nonuniversal amplitude 
factors, universal scaling functions. The scaling expo- 
nents of the correlation length and the crossover variable 
are 



1 



2 - Kl 

1 - K2 



11 

450^ 



2-.r2-45^ + ^(^^)- 



(4.37a) 
(4.37b) 



In the following we will concentrate on a path in the 
phase diagram spanned by the relevant variables r and 

V which approaches the tricritical point r = ?; = as 

V ^ T. Hence, we will neglect the crossover variable 



v/\a\ 



,1-0 



<C 1, where we have set fJ. — 1 for con- 



venience. In this regime we can write the fundamental 
scaling form H4.36|l as 

G^,^({r},r) = \rf^+^^' F^^ahrv}), (4.38) 



with 



,3 = „(d-4 + 5) =i-gE + 0(£^). (4.39a) 



(4.39b) 



The superscript of the scaling functions F^ ^ corre- 
sponds to the sign of t. Note that ^ ~ 1^1" is the 
correlation length. 



2. Scaling behavior of various percolation observables 

Now we will exploit our knowledge about the correla- 
tion functions of the flelds to extract the scaling behavior 
of various observables that play an important role in per- 
colation theory. First we will consider the case that the 
process starts with a single seed at the origin r = 0. Sec- 
ond we will look at the case that the density of the initial 
state is homogeneous. 

Let us start by considering clusters of a given finite 
size S, i.e. clusters with a finite mass S of the debris 
given that the process started with a seed, a single agent, 
at the origin r = 0. In principle, any reasonable ini- 
tial state can be prepared by choosing the appropriate 
seed density po{v). In the Langevin equation (|3.1() this 
general initial condition corresponds to an additional 
source term A"^g(r,t) = po{r)S{\t). At the level of 
the dynamic response functional (|3.7|l such an ini- 
tial state translates into a further additive contribution 
— J d'^x dt q{r,t)h{r,t). Thus, a seed q{r,t) ~ S{r)S{t) is 
represented by the contribution — ri(O,0). At the level 
of the quasistatic Hamiltonian H (|3.11|l such a seed is 
therefore represented by an additive term — s(0). 

Let P{S)dS be the measure for the probability that 
the cluster mass of the debris generated by a seed at the 
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origin is between S and S + dS. In our field theoretic for- 
mulation the probability density P{S) can be expressed 



P{S) ^(^^[J d'^r s(r) - S*) exp [5(0)]^ . (4, 



40) 



For big clusters with S' ^ 1 we can expand the exponen- 
tial to first order (higher orders lead asymptotically only 
to subleading corrections) and obtain 



Pas (5) = ^'^(y d'^J 



s{v)-s \ m 



(4.41) 



for the asymptotic distribution. We will return to Pns{S) 
in a moment. 

The percolation probability Voo is defined as the prob- 
ability for the existence of an infinite cluster generated 
from a single seed. Hence Poo is given by 

/•OC 

V^^l- lim / dS'e-"^P(S') 

= 1- \im^{exp (l{0) ~ c J d'^r s{r)^) . (4.42) 

Via expanding exp[s(0)] we obtain the asymptotic form 
M 



Voo^- lim (S(O)c- 



cM\ 



(4.43) 



where M = J d'^x s{x.). The virtue of this formula is 
that it relates the percolation probability unambiguously 
to an expression accessible by field theory. For actual 
calculations the term exp(— cAf ) has to be incorporated 
into the quasi-static Hamiltonian. This leads to 



nc = n + J d'^rc{r)s{r) 



(4.44) 



instead of the original H. Here, c(x) = c is a source 
conjugate to the field s. Whereas in general (s) = if 
c = by virtue of causality, the limit c +0 leads to a 
nonvanishing order parameter Voo in the spontaneously 
symmetry broken active phase. Having introduced He, 
we can write 

Poo = - lim (S(0))c 

= -Go4(0,r,c^ +0,u,Ai), (4.45) 

where (•••)£ denotes averaging with respect to He- With 
the help of Eq. H4.45|l and the scaling form l|4.38|l we 
readily obtain that 



9{-r)\rf' 



(4.46) 



In order to examine the scaling behavior of P{S) we 
can look at its moments defined by 



/>oo 

(5*^) = / dSS''P{S). 
Jo 



(4.47) 



Using Eq. (|4.41() . our scaling result H4.38|l leads to 

(5^)^ J{d''r)''GkAM,0,r)^\Tf-''^''-^'K (4.48) 

This tells us that Pas scales as 

PUS,t) = Snsir) = S^'^" /{tS''-) , (4.49) 

where ns is the number of clusters of size S per lat- 
tice site. The ns play an important role in percolation 
theory where they are called cluster numbers. The scal- 
ing exponents in the scaling form (|4.48|) reflect the usual 
nomenclature of percolation theory. They are given by 



1 



1 



(4.50a) 
(4.50b) 



These exponents coincide with the corresponding expo- 
nents of conventional isotropic percolation only in mean- 
field theory. It follows from Eq. (|4.48() that the mean 
cluster mass (S) of the finite clusters scales as 



{S)^M{t)=Mo \r\-\ 
with the exponent 

7 = d^.-(/3 + /?')-l + ^e + 0(£2), 



(4.51) 



(4.52) 



Next we consider Green functions restricted to clusters 
of given mass S. In terms of the conventional unrestricted 
averages with respect to H, these restricted Green func- 
tions can be expressed for large S as 



Cf\{r},r) 



(rl)•••s(r^^)5( / rf-^r s(r) - 5)?(0) 



(conn) 



(4.53) 



Equation 



leads to the scaling form 



C^f )({r}, r) = Fw({|rr r}, Irf"-^ S) . 

(4.54) 

With help of these restricted Green functions we can 
write the radius of gyration (mean-square cluster radius) 
of clusters of size S as 



/d''rr2Cp'(r,r) 
2d J d'irC[^\r,T) ' 

Equation (|4.54() then leads to 

with the fractal dimension 

Df=d-^-=^-^,e + Ois^) 



(4.55) 



(4.56) 



(4.57) 
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We conclude Sec. lIVI bv considering the scaling behav- 
ior of the statistics of the debris if the initial state is 
prepared with a homogeneous seed density po. As dis- 
cussed above, such an initial state translates at the level 
of the quasistatic Hamiltonian H (|3.11|l into a further 
additive contribution — po J d'^r's{r). Our general scal- 
ing form (|4.38|l implies that the correlation functions of 
the densities s(r) scale in case of the homogeneous initial 
condition like 

G^,^({r}, r, po) = E § / (^'^' G^^^^,{{r}, r) 

k=0 ■ 

= kr ^jv({kl po) ■ 

(4.58) 

It is obvious that the initial seed density plays the role 
of an ordering field. Hence, the Green functions do not 
show critical singularities as long as po is finite. For the 
homogeneous initial condition the appropriate order pa- 
rameter is given by the density of the debris 



P= (s(i")>po Gi,o(r,r,po) . 
Equation H4.58|) tells us that 



P = 



\rff^i\ 



rr Po)- 



(4.59) 



(4.60) 



In the non-percolating phase (r > 0) the order parameter 
p is linear in po for small seed density with a susceptibility 
coefficient that diverges as t ^ 0, 



p{t > 0,po) ''po; 



(4.61) 



with the susceptibility exponent 7 given in Eq. H4.52|l . 
At criticality (t — 0) the order parameter p goes to zero 
for Po ^ as 



p(t = 0,po) ^ Po 



1/5 



with the exponent 



= 1 



(4.62) 



(4.63) 



Finally, in the percolating phase (r < 0) the order pa- 
rameter becomes independent of the initial seed density 
in the limit po ^ and goes to zero with r as 



p(r < 0) 



(4.64) 



Equations (|4.46|) and (|4.64|) show explicitly that the 
two order parameters, namely the density p of the de- 
bris and the percolation probability Voo^ have different 
exponents [3 and /?'. This stands in contrast to ordinary 
isotropic as well as directed percolation. As mentioned 
earlier, ordinary isotropic and directed percolation are 
special in the sense that they posses an asymptotic sym- 
metry that leads to the equality of the respective expo- 
nents (3 and (3' . 



E. Scaling properties in 5 dimensions — logarithmic 
corrections 

Here we will study TIP directly in 5 dimensions where 
fluctuations induce logarithmic corrections to the lead- 
ing mean-field terms rather than anomalous exponents. 
First, we will establish the general form of the logarith- 
mic corrections. Second, we will analyze how these cor- 
rections, to leading order, affect the observables studied 
in Sec. UVDl 



1. General form of the logarithmic corrections 

Now we will solve the characteristics directly for rf = 5. 
The Wilson functions stated in Eqs. (|07|) and 
are central ingredients of the characteristics. For cc- 
conomic reasons, we will use in the following an ab- 
breviated notion for the Wilson functions of the type 
J{u) — fiu + /2U^ + • • • . For instance, we will write 
Eq. (|4.28a|l as /3„(u) = Piu + l32u'^ + 0{u'^) and likewise 
for the other Wilson functions. Since we are interested 
in the tricritical point, we set v = throughout. 

Solving the characteristic for u, Eq. H4.30|l . for e = 
yields readily 



e = e{w) = io exp 



J- + 0{lnw) 

P2W 



(4.65) 



where io is an integration constant and where we abbre- 
viated 

w = u . (4.66) 

To leading order we obtain from Eq. H4.65|l asymptoti- 
cally for ^ < 1 



1 



45|ln^| 



< 1. 



(4.67) 



The remaining characteristics, Eq. H4.31|l are solved via 
exploiting id/d£ = Pd/dw with the asymptotic result 



Q(u;) ~ Qow«l^^^ 



(4.68) 



where Qo symbolizes a non-universal integration con- 
stant. Having solved the characteristics, we obtain 



G^^({r},r) = [^,£Xiwy/Y[WX{wy/Y 



X F 



{Pir} 



t{w) 



(4.69) 



as a formal scaling solution for the Green functions. 

In the following we have to be careful about the ex- 
plicit dependence of the scaling functions Fj^ on w 
even though we are only interested in the leading log- 
arithmic corrections. This intricacy comes from the fact 
that rather mean-field than Gaussian theory applies at 
the upper critical dimension. Hence we must carefully 
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distinguish between the 2 roles of g, viz. its role as a 
dangerous irrelevant variable which scales the fields and 
the parameters in the correlation and response functions, 
and, its role (in its dimensionless form u) as the loop ex- 
pansion parameter. Only the latter role can be safely 
neglected for the leading logarithmic corrections. As we 
shall see in the following, the dependence of the functions 

Fn.n ^ if gi'^™ by 



F'^~{{r},T,v/V^)+0{V^) . (4.70) 



To derive this result one can consider the loop-scaling of 
the generating functionals of the vertex and Green func- 
tions as it was done in Ref. js^ for studying logarithmic 
corrections in DP. Here we use a more direct method to 
determine the dependence of the tree diagrams on the 
dangerous variable g. Consider an arbitrary connected 
diagram with L loops, P propagators, N external s-legs, 
N external s-legs, Vg vertices of type g, Vf vertices of type 
/, and Vi noise vertices proportional to 1. This generic 
diagram contributing to the Green function Gjy jy satis- 
fies the three topological conditions 

L = P-Vg-Vf-Vi-N-N+l, (4.71a) 
P = 3Vg + 2Vf + Vi+N, (4.71b) 

P = Vg + Vf + 2V1 + N- (4.71c) 

Eliminating P we arrive at two equations for the number 
of vertices, namely 

Vi=N-l + L, (4.72a) 
2Vg + Vf = N -1 + L. (4.72b) 

Switching from / to the scaled variable /' = f/y/g, we 

find that the diagram scales with g as with 
the factor g^^^ being determined by the loop-order of the 
diagram. Upon renormalization this reasoning leads for 
L = to the scaling form (|4.70|) . 

Knowing Eq. H4.70|l we can write down the scaling form 
of the Green functions with its leading logarithmic cor- 
rections. From the general solution of the RGE for the 
Green functions (|4.69|) in conjunction with the asymp- 
totic result H4.67|l we get 

Gj, ^({r} , r) = ^^+2* |i^^^|77V/75-22Ar/225+(l-7V)/2 



where f-p^ is non zero only for r < 0. Now we fix the ar- 
bitrary but small flow parameter £ so that effectively 
acquires a finite value in the scaling limit. 



|ln^| 



-11/225 , 



1. 



Hence, we choose asymptotically 



|ln|r| 



1-11/225 



Collecting, we obtain that 

rM^e(-r)\T\\\n\T\f''. 



(4.75) 



(4.76) 



(4.77) 



Next we consider the probability P{S,t). Using the 
definition H4.48|l and the general result H4.73f) we obtain 

X /},(r2|ln^r"/''^T,£4|ln^|22/225 ^) (4 78) 

Being interested primarily in P{S, r) as a function near 
criticality r « and not in the animal limit S oo with 
T > 0, we hold the second argument finite. Hence we 
choose 



^4^^-l(ln^)-22/225^ 

With these settings we find 

P{S, t) = S-^'^ (In 5)2/45 fp (r5i/2) . 



(4.79) 



(4.80) 



The scaling function fp (x) is found with help of a 
simple mean field calculation to be proportional to 
exp(— const x^) as long as x is finite. For x ^ oo it crosses 
over to the animal limit with its own independent scaling 
behavior. Note that there is, to the order we are work- 
ing, no logarithmic correction associated with the S'^/2 
in the argument of the scaling function. As a corollary 
of Eq. (|4.80(l , we obtain the critical behavior of the mean 
cluster mass, 



M{t) = {S) 



-1 n |2/45 



(4.81) 



Now we turn to clusters of a given size S. The scal- 
ing behavior of the Green functions restricted to these 
clusters takes on the form 

4^'({r},T)=^^+6|ln^|(^^-22^)/22^ 



X F^,^({M,^^'|ln^r"/'''r). (4.73) ^ F^({M llnff"/^^^ r,^^ |i^^|22/225 (482^ 



2. Logarithmic corrections to the percolation ohservahles 

As above we will first consider the case that the process 
emanates from a single local seed at the origin. Exploit- 
ing Eq. (|4.45() . we find that the percolation probability 
has the scaling form 

Voo=f\\ni\'"'fv^{r^\ln£\-^^'^^''T), (4.74) 



in 5 dimensions. From this scaling it is straightforward 
to extract the radius of gyration of clusters of size S via 
definition H4.55|l . Choosing t per Eq. (|4.79() leads to 



Rl^S^/^ (ln5)"/225/j,(r5i/2). 



(4.83) 



In mean-field theory, as a straightforward calculation 
shows, the scaling function fB.{x) is identical to a non- 
universal constant. 



In the remainder of Sec. lIVEl we will consider the ho- 
mogeneous initial condition. The A'^-point density corre- 
lation functions obey in d = 5 the scaling form 

(4.84) 

Specifying Eq. H4.84|l to = 1 and iV = 0, we obtain the 
scaling behavior of the mean density of the debris, 

X Fi,o(r2|ln£r"/22^r,r3|ln^r«i/i^%o). (4.85) 

For r 7^ and small seed density po it is appropriate to 
choose t according to Eq. (|4.75(l . This provides us with 

-p{r^p,) = \r\'" Mr\f"'ff{p, |lnM|-/^). 

(4.86) 

The scaling function jf{x) behaves like f^{x) ~ x for 
small X, and hence 

p(r>0,po)~Por-MlnT|'/"' (4.87) 

in this regime. At t = 0, the mean density is a function 
of po only. To obtain the logarithmic corrections for this 
case, we set 

l^p,^/^ \\np,\-'''''r (4.88) 
From Eq. H4.85|l in conjunction with Eq. (|4.88|) we obtain 

-p[T^Q,p,)^pl"\\np,t"'' . (4.89) 

Lastly, the order parameter should be independent of the 
seed density in the percolating phase for po ^ 0, i.e., 
the scaling function fp of Eq. H4.86|l should approach a 
constant in this limit. Thus, we get 

P(t < 0) - |r|'/' |ln|r||"/^^ (4.90) 

in this regime. 

V. DYNAMIC SCALING PROPERTIES 

In this section we will study dynamic scaling proper- 
ties of TdlP. To find out these properties we have to 
renormalize the dynamic response functional (|3.7|) . In 
comparison to the quasi-static Hamiltonian H3.11|l , has 
one additional parameter, namely the kinetic coefficient 
A. To determine the renormalization of A, we have to 
calculate one of the dynamic vertex functions in full time 
or frequency dependence. A dynamic RGE then leads 
to general scaling form for the dynamic Green functions. 
This scaling form allows us to deduce the dynamic scaling 
behavior of various percolation observables. 
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(a) (b) 

FIG. 10: Dynamic two-loop diagrams which provide us with 
the renormalization factor of the kinetic coefficient A. 

A. Diagrammatics 

In order to determine the renormalization of A, we 
could calculate the frequency dependent part of any of 
the superficially divergent vertex functions. For conve- 
nience, we choose to work with Fi^i. Two dynamic two- 
loop diagrams for the self-energy can be constructed from 
the diagrammatic elements in Fig.|21 These two diagrams 
are shown in Fig. 1101 After Fourier transformation and 
some rearrangements, we get from the diagrams the con- 
tribution 

Cnt) + CI3') = iu:gT-' ('_| + | + ^/,22 + ^/,,3') (5.1) 



to Fi i, where we have not displayed various terms (those 
not linear in the frequency uj) for notational simplicity. 
B and C stand for the integrals 




(5.3) 

These integrals can be calculated for example by using 
Schwinger parametrization, cf. the appendix. The e ex- 
pansion results for B and C read, 

i?=^(3^/3-.), C = (5.4) 

The /...-contributions in Eq. H5.1() cancel and we finally 
get 

m) + m) = — f — + (5.5) 

B. Renormalization and renormalization group 
equation 

The dynamic theory requires an additional renormal- 
ization factor, say Z', in comparison to the quasi-static 
theory due to the existence of e poles in the frequency 
dependent terms. Taking into account that the dynamic 
fields n and n are related to the quasi-static fields s and s 
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via Eq. (|3.10() , we introduce Z' consistent with the renor- 
malization scheme in Eqs. (|4.16|l and (|4.17l) by letting 



- Z'l/2, 



(5.6) 

Combining Eq. (|5.5|) with the zero-loop part of the vertex 
function Fi i we get by applying our renormalizations 



(Z'Z)^/^-ii(l + l^)+0(.^) 



(5.7) 

To keep the formula H5.7|l simple, we have once more 
dropped all terms which are not linear in lo. From 
Eq. I|5.7|l we can directly read off (Z'Z)^/^ to order u. 
Taking into account Eq. H4.20|l . this yields 



^' = i + -f^-?^Vo(.^) 

e \ TT 15/ 
The corresponding Wilson function reads 



/58_4-V3 

' V15 TT 



(5.8) 



(5.9) 



Now we have all the information required to calculate the 
Gell-Mann-Low function 

C = A^a^lolnA-l(7'-7) (5.10) 

for the kinetic coefficient. 

By proceeding analogous to the quasi-static case we 
obtain the dynamic RGE 



I?^ + i(iV7' + iV7) 



G^^({r,i},T)=0. (5.11) 



Here, the RG differential operator I?^ is given by 

(5.12) 

C. Scaling properties 

1. General scaling form 

Upon using the method of characteristics we obtain 
the dynamic scaling form 

G^,^({r,0,r)=/-^^('^-2)^+^^ 

X Fj, ^({^/ir, e^Xfih}, fi-^a/e'/"-, fi-^/e't'/n , (5.13) 



where 



S^j,^[d-A+^ + z)N+[2 + ^)N 



= (- + z)n+^-N. (5.14a) 

„ „ A 8 4- \/3\ 2e ^, 2n 

^-2 + ^ 15+^ 45+^(^) 



2 - 0.0558 £-fO(e2). 



(5.14b) 



The dynamic exponent z is identical to the fractal di- 
mension Dmm of the minimal or chemical path, Dmin = 
z 33]. Near the tricritical point [v = 0), we get from 
Eq. (|5.13|) that the response and correlation functions of 
the agent at time t obey the scaling form 

^({r}, t, r) = ^({rA^/^, rtV- ) , 

(5.15) 



' N,N 

where 



VZ \ 6 TT / 45 



^"^^f + 0(£^), (5.16a) 



= VZ = 1 



' +0{e'). (5.16b) 



TT / 45 



2. Dynamic scaling behavior of various percolation 
observables 



First we will consider the spreading of the agent ema- 
nating from a localized seed at r = and t — 0. Later 
on, we will turn to the case that the initial state at time 
t = is prepared with a homogeneous initial density po- 

The survival probability V(t,T) that a cluster grown 
from a single seed is still active at time t can be derived 
from the field theoretic correlation functions by using [s^ 



P{t) = - lim (e-'=^(°)n(-t)) , 



(5.17) 



where now A^(0) — J d'^x n{x,0). By proceeding analo- 
gous to the static case, i.e., by incorporating the term 
exp(— fcAf(O)) into the dynamic functional, one obtains 



V{t,T) = - lim (S(-i))fe 

A;— >oo 



-G'o,i(— ^, r^k — cxD, u; A, fi) . 



(5.18) 



where {■ ■ ■)k denotes averaging with respect to the new 
dynamic functional J7fc that has absorbed the source term 
featuring k. Equation H5.15|l then implies that the sur- 
vival probability obeys the scaling form 



(5.19) 



In the percolating phase (r < 0) the survival probability 
tends to the percolation probability for large t, 'P{t, r) 

Vooij) ~ |r|^ . Hence, the universal scaling function 
fv{x) behaves as for large values of x. 

The mean density of the agent at time t grown from a 
single seed follows from Eq. H5.15|l with = = 1 as 

p{r, t, r) = ^-l-(^+'5')/'^V(^A^^^ Ti^^"^) • (5.20) 

Knowing the density of the agent, we have immediately 
access to the number of infected or growth sites 



(5.21) 
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which can be viewed as the average size of the epidemic 
at time t. Equation H5.20|l impHes that 



where 



r/. - — - 1 
vz 



1 4-V3\ e 



(5.22) 



TT / 45 



+ 0(e^). (5.23) 



In the non-percolating phase (r > 0), the integral 
/p dtM(t,T) is proportional to the mean mass {S) of 
the static clusters of the debris, Eq. (|4.51(l . The mean 
mass of the debris at time t, on the other hand, is given 

by 



(5.24) 



(5.25) 



M{t,T) = / dt'JV{t',T). 
Jo 

For M we find the scaling form 
with 

7?, = 1 + 7], = ^ = ^ . (5.26) 

VZ 

The scaling functions ju(x) and /ai(x) are regular for 
small X. For a; — *■ oo we learn from limt^oo -/^(^j ''') — 
M.{t) = (S) in conjunction with Eq. H4.51|l that 



Ma - fj^iy) - exp(-const?;). 



(5.27a) 
(5.27b) 



It follows that the number of growth sites behaves for 
r > asymptotically as 



Af{t,T) - T'^^-'^exp(-constT'-'^i). 



(5.28) 



Knowing the density of the agent, we are in the posi- 
tion to calculate the mean square distance Tl^{t, r) of the 
infected individuals from the original seed by using 



1 



ArV(r,i,T). (5.29) 



2dAf{t,T) 
We obtain the scaling behavior 

7^2(^,r) -^^=/7^(T^l/''0, 

where 



(5.30) 



z \ 15 TT / 45 

Mendes et al. proposed a generalized hyperscaling 
relation which translates in our case to 



2[l + ^)5, + 2f],^dzs 



(5.32) 



From Eqs. (|4.51() . \b.l9\ . H5.25|l . and H5.29|l we we confirm 
that the spreading exponents indeed fulfill this relation. 
This is not a surprise because the hyperscaling relation 
(|5.32l) is based only on the general scaling form (|5.15|) . 

Finally, we consider the scaling behavior of the time de- 
pendent mean density of the agent p{t, r, po) — ("-(r, t)) pa 
for i > if the initial state at time i = is prepared with 
a homogeneous initial density po- As mentioned ear- 
lier, this initial condition corresponds to a source term 
q(r,t) = poS{t) in the Langevin equation (|3.1|l . This 
source term translates into a further additive contribu- 
tion —po J d'^rnir, 0) to the dynamic functional H3.7|l . In 
a theory like ours, where the perturbation expansion is 
based only on causal propagators and where no correla- 
tors appear, no initial time UV-infinities are generated. 
Therefore, no independent short time scaling behavior 
[27LI3H arises and n(r, 0) scales as n(r, t). Thus, we find, 
analogous to Eq. I|4.58|l . that the dependence of the cor- 
relation functions on po can be expressed as 



Cj,i{r},t,r,po)^t-^'+^/'^^^'' 



(5.33) 



In particular, we obtain for the mean density of the agent 

p{t, T, po) = ^-l-'5/^Vp(r^^/'^^ Poi'"''-''')/'^") . (5.34) 

At criticality (t = 0) it follows from this equation that 
the agent density first increases in the universal initial 
time regime, 

p{t,po)- Pnt''' . (5.35) 
Then, after some crossover time, it decreases, 

p(t,Po)^i-'^^/'^% (5.36) 
with the critical exponent 



1 + /107 4-V3 \ £ 



uz 2 \ 3 



TT / 90 



+ 0(^2) . (5.37) 



The time dependence of the order parameter, the den- 
sity of the converted individuals 

P(i, T, Po)=/ dt' p{t',T,po), (5.38) 
Jo 

follows from Eq. (|5.34(l . 

p{t, T, po) = t-''/''^fp{Tt'/-'\pot^'"'-'^^/-'^) . (5.39) 

This scaling form goes exponentially to the time indepen- 
dent scaling form (|4.60|) in the large time limit. Equa- 
tion ()5.39|) implies 



p(i,0,po) Pot\ 



(5.40) 



for the initial time order parameter scaling at criticality. 
The scaling index appearing here is 



J_ 
vz 



(5.41) 
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D. Logarithmic corrections in 5 dimensions 

Here we are going to investigate how logarithmic cor- 
rections influence the dynamic scahng behavior in = 5. 
We win briefly explain how the general considerations on 
logarithmic corrections in the quasi-static theory given 
in Sec. IIVETI have to be augmented and modifled in the 
dynamic theory. Then we will derive the logarithmic cor- 
rections for all of the dynamic observables studied above. 

1. General form of the logarithmic corrections 

Compared to the quasi-static theory, X' takes on the 
role of X and there is an additional flowing variable, viz. 
A. Both X' and A have characteristic equations of the 
form given in Eq. H4.31|l and hence the flow of these vari- 
ables in d = 5 is described by Eq. (|4.68(l . The dynamic 
scaling form H5.13|l becomes 

G^^^({r,t},r) = i^^lf^'+'^X'iwf/'Xiwf^' 

X ^w,7v( {/^^r, {fie)^X{w)t} , ^) • (5.42) 

in 5 dimensions. Since we are interested solely in the 
leading logarithmic corrections, it is sufficient for our pur- 
poses to account for the explicit dependence of the scal- 
ing functions on to 0-Loop order, see Eq. 14.70|l . 
Using the asymptotic result H4.67|l we obtain 

g3N+2N n^^|77V/75+(a-22/225)W+(l-7V)/2 

xFj,,^{{£r,f\\ni\''t}J-'\ln£\''/'^'T), (5.43) 
where a is the abbreviation of 

As it was in Sec. IV C 21 our emphasize here is on the 
time dependence of various measurable quantities. Hence 
we fix the arbitrary flow parameter £ as long as A |t| t < 1 
by setting 

^2 ^ (ln^)-^ (5.45) 
In the contrary case A |r| f ^ 1 we must use the choice 

(Ezni). 

2. Logarithmic corrections to dynamic percolation 
observables 

As above we will first consider the initial condition that 
the process starts from a single local seed at the space- 
wise and time-wise origin. The first quantity that we 
are going to consider is the survival probability V{t,T). 



Utilizing Eq. H5.18|l in conjunction with Eq. (|5.42l) we 

obtain 

V{t, t) = e lln^l^/^^ Yv{£' |ln£|" t, jln^p"/'"^ r). 

(5.46) 

The choice (|5.45ll then leads to 

V{t, t) - {\ntr^ friTt{\ntr-), (5.47) 

with the exponents 

7 1 /II 4- x/3\ 
a, = --a = -(- — ), (5.48a) 

a.=a-ii^lf-Ui^). (5.48b) 
225 45 V 3 TT / ^ ' 

Hence, we have at criticality and near criticality with 
A|r|t < 1 

r{t,T^O) ^t-'^{\nty^. (5.49) 

The asymptotic behavior for large times below and above 
criticality, respectively, is given by 

V{t,T> 0) - T llnrl^^^^ exp(-constir |lnr|"^) 

(5.50a) 

and 

rit,T <o)-rooi-\T\) 

- |t| |ln|T||'^^^ exp(-consti|T| llnlrH"^). 

(5.50b) 

Next we look at the mean density of the active par- 
ticles at time t. Upon specializing the general scaling 
form to iV = TV = 1 we find 

p(r,t,r) = £^ lln^l"^ fpi^rj^ \\ii£\U,£-^ |ln£r"/''^r) 

(5.51) 

with 

From the mean density we obtain the mean number of 
agents (|5.21l) at time t by integrating over r, and choosing 
01^ for T « 

Mit,T) = \\ntr f^{Tt{lntr^) . (5.53) 

For the mean mass of the debris at time t as defined in 
Eq. H5.24|l we get 

X(i,r) =t|lntr^/M(ri(lnt)°-). (5.54) 

Ai (t, r) crosses over for r > and t ^ oo to a function 
which approaches (t) — (S), Eq. (|4.81|) . exponentially. 
The mean square distance TZ^ {t, t) of the agents from the 
origin is found to behave as 

n^it^T) ^tilntr fn{Tt(lntr-). (5.55) 
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Finally, we will consider the homogeneous initial con- 
dition. The analog of Eq. H5.33|l in 5 dimensions with the 
flow parameter £ still arbitrary reads 

^ /r 1 , N /j3N \^ /,|(a-22/225)JV+l/2 

X ^^Ar({^r,£2|i^^|a^|^^-2|i^^|-ll/225 ^^^-3|i^^|-61/150^^^ 

(5.56) 

where we have used Eq. H5.43|l . From Eq. (|5.56() we read- 
ily obtain the critical behavior of the mean density of the 
agents by setting = 1 and fixing £ for not too large t 
via Eq. 15.45|l . 

pit,T,po) = (Int)-p fp{rt {lntr-,pot'/' (Int)-""), 

(5.57) 

where 

450 2 90 V 3 TT /' ^ ' 

a^ = ^~^ = 1(^-1^). (5.58b) 
150 2 30 V 3 IT J ^ ' 

At criticality, the scaling function is expected to behave 
as fp{0,y) ~ y for y ^ 1. Thus, the agent density in- 
creases initially, 

p(i,0,po)-po(lnt)"^. (5.59) 

After some crossover time it decreases like 

pit,0,po)^t-^^^ilntr^ . (5.60) 

The mean density of the debris at time t can be extracted 
without much effort by integrating over Eq. H5.57|l . This 
yields 

p{t, T, po) = (ini)a, f^^^^ {\ntr^,po ^3/2 (Int)-'^") 

(5.61) 

if t is not too large. In the case po t^l"^ ^ oo we have to 
use 

^3^po|lnpor''/''°, (5.62) 
which implies especially at criticality 

p(t,r = 0,po)=p;,/'|lnpo|'/'' TX^pT llnpoP''^"/'). 

(5.63) 

The scaling functio n /^( x) approaches exponentially the 
stationary density H4.89|l . 

VI. CONCLUDING REMARKS 

In summary, we have generalized the usual GEP by 
introducing a further state in the live of the individu- 
als governed by the process. Our GGEP has a multi- 
dimensional phase diagram featuring two surfaces sepa- 
rating endemic and pandemic behavior of the epidemic. 
One of the surfaces is a surface of first order phase transi- 
tions whereas the other surface consists of critical points 



representing second order transitions. The two surfaces 
meet at a line of tricritical points. 

The second order phase transitions belong to the uni- 
versality class of dynamic isotropic percolation (dIP). In 
the vicinity of these transitions, the asymptotic time limit 
of the GGEP is governed by the critical exponents of 
usual percolation. The debris left behind by the process 
forms isotropic percolation clusters. 

Mainly, we were interested in the tricritical behavior of 
the GGEP. We set up a field theoretic minimal model in 
the form of a dynamic response functional that allowed 
us to study in detail the static and the dynamic scaling 
behavior of the universality class of tricritical dynamic 
isotropic percolation (TdIP). In particular we calculated 
the scaling exponents for various quantities that play an 
important role in percolation theory. As expected, these 
exponents are different from the exponents pertaining to 
dip. For example, we computed the exponents /3 and /3' 
respectively describing the two different order parame- 
ters, viz. the density of the debris and the percolation 
probability. Whereas /3 and /3' are identical in dIP, they 
are different in TdlP. Although TdIP is described by scal- 
ing exponents different from those of usual percolation, 
we learned that its spreading as well as the statistics of 
its clusters behave in many ways like conventional dy- 
namic percolation. For example, TdIP has meaningful 
cluster numbers, fractal dimensions etc. Thus, we pro- 
pose to refer to the static properties of the TdIP as tri- 
critical isotropic percolation (TIP). 

The surface of first order transitions is characterized 
by a compact cluster growth, i.e., the fractal dimension 
of the clusters is identical to the their embedding dimen- 
sion. We hope that our findings trigger numerical work 
with the aim to verify the predicted first order perco- 
lation transitions. A promising strategy that avoids a 
cumbersome detection of jumps in the order parameters 
might be to measure directly the fractal dimension of 
clusters near the first order surface. 
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APPENDIX: CALCULATION OF THE 
PARAMETER INTEGRAL I{a,b;c) 

In this appendix we sketch our calculation of the pa- 
rameter integral I{a,b;c) defined in Eq. H4.6|l . Most of 
the integrals that have to be performed in calculating the 
two-loop diagrams can be derived from /(a, 5; c) simply 
by taking derivatives with respect to the parameters a, 
b, and c. The integrals B and C, which occur in the dy- 
namic calculation and cannot be extracted from /(a, b; c) 
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by taking derivatives, can be calculated by similar means 
as /(a, b; c). 

In the following we use the so-called Schwinger 
parametrization which is based on the identity 

1 I 

— = —— / dsT-^expf-As) , Ren>0. (A.l) 

A" r(n) Jo J' \ j 

In this parametrization, Eq. H4.6|l takes on the form 

//•oc 
/ dsids2ds3 S3 exp { — si [q? + a] 

qi,q2 

- S2 hi + b]- si [(qi + q2)2 + c] } . (A.2) 

A completion of squares in the momenta renders the mo- 
mentum integrations straightforward. We obtain 



/(a, 6; c) 



1 



i^T^r Jo 

exp(— sio — S26 — S3C) 

[SlSs + S2S3 + SiS2]'*/2 ■ 



(A.3) 



Changing integration variables, si — > tx, S2 — *■ ^(1 — x) 
and S3 tz, and carrying out the integration gives 



/(a, 6; c) 



r(4 - d) 
{^^Y JO 



dz / dx 



X [ax -|- b{l ~ x) + cz 



+ 



2 

(A.4) 



The remaining integrations can be simplified by letting 
z {z~^ — l)x{l — x). After this step, which leads to 



I(a, b: c) = — ^- — -7^ / dz [ dx 
(4^)'^ Jo Jo 



X [axz + 6(1 — x)z + cx{l — x){l — z] 



(A.5) 



one sees easily that the remaining integrations are finite 
at the upper critical dimension. Hence, they can be con- 
veniently evaluated directly at d = 5. An e expansion of 
the gamma function, 



r(4-d) = -£ii±i/^ + 0(e"), (A.6) 



finally leads to the result for I {a, b; c) stated in Eq. (|4.7|l . 
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